3D Euler about a 2D Symmetry Plane 



Miguel D. Bustamantc 1 ' * and Robert M. Kerr 1 
1 Mathematics Institute, University of Warwick, Coventry CV4 7AL, United Kingdom 

Initial results from new calculations of interacting anti-parallel Euler vortices are presented with 
the objective of understanding the origins of singular scaling presented by Kerr (1993) and the lack 
thereof by Hou and Li (2006). Core profiles designed to reproduce the two results are presented, new 
more robust analysis is proposed, and new criteria for when calculations should be terminated are 
introduced and compared with classical resolution studies and spectral convergence tests. Most of 
the analysis is on a 512 x 128 x 2048 mesh, with new analysis on a just completed 1024 x 256 x 2048 
used to confirm trends. One might hypothesize that there is a finite-time singularity with enstrophy 
growth like __ ~ (T c — t)~ ln and vorticity growth like ||<a>||o_ ~ (T c — i)~ 7 . The new analysis would 
then support 7n w 1/2 and 7 > 1. These represent modifications of the conclusions of Kerr (1993). 
Issues that might arise at higher resolution are discussed. 
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I. INTRODUCTION 

One definition of solving Eulcr's three-dimension in- 
compressible equations [1] is determining whether or not 
they dynamically generate a finite-time singularity if the 
initial conditions are smooth, in a bounded domain and 
have finite energy. The primary analytic constraint that 
must be satisfied [2] is: 

f IMUdt-00 (1) 
Jo 

where ||u>||oo is the maximum of vorticity over all space. 
To date, Kerr (1993) [3] remains the only fully three- 
dimensional simulation of Eulcr's equations with evi- 
dence for a singularity consistent with this and related 
constraints [4]. Growth of the enstrophy production and 
stretching along the vorticity, plus collapse of positions, 
supported this claim [3] . Additional weaker evidence re- 
lated to blow-up in velocity and collapsing scaling func- 
tions was presented later [5]. 

There is only weak numerical evidence supporting 
these claims [6, 7]. In a recent paper, as described in 
one of the invited talks of this symposium, Hou and Li 
(2006) [8] found evidence that the above scenario failed 
at late times. 

This contribution will first comment on four issues 
raised at the symposium, then present preliminary new 
results. The four issues are: 

• How should spurious high-wavenumber energy in 
spectral methods be suppressed? 

• What criteria should be used to determine when 
numerical errors are substantial? 

• What effect do the initial conditions have on singu- 
lar trends? A cleaner initial condition is proposed. 
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• We introduce a new approach for determining 
whether there is singular behavior of the primary 
properties and the associated scaling. This is ap- 
plied to both new and old data. 

All calculations will be in the following domain: L x x 
L y x L z = 4-7T x 4-7T x 27r with free-slip symmetries in 
y and z and periodic in x with up to n x x n y x n z = 
1024 x 256 x 2048 mesh points. Using these symmetries 
only one-half of one of the anti-parallel vortices needs to 
be simulated. 

The "symmetry" plane will be defined free-slip 
symmetry through the maximum perturbation of the ini- 
tial vortices and the "dividing" plane will be defined as 
the xy free-slip symmetry between the vortices. 



II. HOW SHOULD SPURIOUS 
HIGH-WAVENUMBER ENERGY IN SPECTRAL 
METHODS BE SUPPRESSED? 

A generic difficulty in applying spectral methods to lo- 
calized physical space phenomena is the accumulation of 
spurious high-wavenumber energy that leads to numeri- 
cal errors. 

What is the best approach for eliminating these spuri- 
ous modes? We have compared the old-fashioned 2/3rds 
dcaliasing versus the recently proposed 36th-power hy- 
perviscous filter [8, 11]. Detailed tests to be described in 
a later paper show that the latter is better in the sense 
that for several quantities, such as the peak vorticity, 
lower resolution calculations follow the high resolution 
cases longer. But a combination of the two approaches 
works even better, and that is what is used here. 

Still, caution is required for any of these approaches 
as the hyperviscosity can dissipate small structures such 
as the anomalous negative vorticity in the squared-off 
profile below. Surprisingly the 36th-order hyperviscosity 
does not appear to produce the ghost vortices that are a 
known artifact of lower-order schemes. 
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III. WHAT CRITERIA SHOULD BE USED TO 
DETERMINE WHEN NUMERICAL ERRORS 
ARE SUBSTANTIAL? 

There are traditionally two approaches to this prob- 
lem, one emphasizing local quantities such as ||k;||oo , and 
the other emphasizing global quantities such as the mean 
square vorticity or enstrophy. We use both. 

1. Local quantities and resolution 

To determine local resolution it is important to check 
the convergence of local quantities such as: 

• The maximum of vorticity ||o;||oo • The location of 
1 1 a; | loo will be defined as x^. 

• The local stretching of vorticity 

a = u>ieijU>j (2) 
where lj = u>/|w| and = ^(uj.j + Ujj). 

Following earlier work [3, 8], we use the criteria that 
Xrx, cannot be closer than 6 mesh points from the dividing 
plane. 

2. Integral quantities 

Examples of integral quantities we could monitor are: 
energy, circulation (which are in principle conserved) , en- 
strophy and hclicity (which are in principle changing) . 

i) Energy is robustly conserved by spectral methods 
even when under-resolved and therefore is not a 
useful test. Convergence of the energy spectrum 
[8] is only a partial test because it neglects phase 
errors. 

ii) Circulation in the upper half of the symmetry plane 
(i.e., the z > half of the x,z-planc, which is 
perpendicular to the primary direction of vortic- 
ity y) is conserved. Circulation in the equivalent 
half of the dividing plane is also conserved. In all 
of the initial conditions considered here, it is ini- 
tially zero and ideally should remain so. There- 
fore, the circulations of the symmetry and divid- 
ing planes, a y = J z>0 u> y (x, 0, z, t)dx dz and a z — 
I y >o L0 ^ x -> V' 0' t)d x dy were monitored. 

We have found that serious depletion of a y is con- 
trolled by n z and the time this begins is indepen- 
dent of the high-wavenumber filter. Once n z is set, 
by convergence of ||u>||oo , we find that there is good 
convergence if n x — n z /2 and n y = n z /A. A later 
paper will provide more details on these conver- 
gence tests. We will violate the condition on n y at 
late times due to current memory restrictions. 



Without the circulation test, it is difficult to draw 
conclusions about the late times in Hou and Li 
(2006) [8] where they claim to see divergence from 
the scaling of Kerr (1993) [3]. 

iii) Enstrophy f2 grows in time, so one test is to check 
how it is balanced by its production Q p , which we 
determine directly. The enstrophy and its produc- 
tion are 

fl = J dVuj 2 , fl p — 2 J dVuJieijLjj . (3) 

iv) Helicity grows within the quadrant simulated (not 
over the full anti-parallel geometry), but its pro- 
duction is determined by pressure which has not 
been calculated. 



IV. WHAT IS THE EFFECT OF THE INITIAL 
CONDITIONS ON THE POTENTIALLY 
SINGULAR BEHAVIOR? 

A. Earlier descriptions 

Because ambiguities in the earlier description of the 
initial condition of [3] led to differences in the initial con- 
dition of Hou and Li (2006) [8] , the community needs a 
clear description of a reproducible, clean initial condition 
that yields the trends of Kerr (1993) [3]. Ideally, we want 
an initial condition whose vorticity is purely positive in 
the upper half of the symmetry plane, which following 
Kelvin's theorem will remain positive for all subsequent 
times. These steps were used [3] to massage the vortex 
profile in order to achieve this: 

1) The first step in creating the initial profile of the 
vorticity core is to use an explicit function where 
the value and all derivatives went smoothly to zero 
at a given radius. See references in [3] for earlier 
work that had used a similar profile. To this, a 
localized perturbation in its position in x was given 
[9]. 

2) The second step is to remove high-wavenumber 
noise by applying a symmetric high-wavenumber 
filter of the form: exp(-a(fc^ + k y + kl) 2 ). Kerr 
(1992) [10] showed the undesirable side-effects if 
this is not done. However, it has become apparent 
that the high-wavenumber filter is not sufficient. 

B. Effect of a negative region 

The upper frames in Fig. 1 come from a reproduc- 
tion of the squared-off profile of Hou and Li (2006) [8] 
which follows the procedure above with the exception of 
using a different high-wavenumber filter. Note the neg- 
ative region in the lee of the primary vortex and how 
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FIG. 1: u)y in the symmetry plane from two initial conditions for t = and for early times with roughly the same growth in 
max(ojj / ) in the symmetry plane. The first squared profile is nearly the same as used by Hou and Li (2006) [8] using step 1) with 
max(ojj / ) = 0.49 in the symmetry plane (over all space ||w||oo = max(|w|) = .67) and for step 2) a squared-off high-wavenumber 
filter (exp(— a(k* + k* + k*)). Note the large negative vorticity in the lee (right) of the primary vortex as in Hou and Li (2006) 
[8] (Fig. 2) and how this is entrained underneath the primary vortex at t = 6, whereupon the hyperviscous filter will dissipate 
it. In the new profile step 3 is included: adding positive ui y {z). In this case max(w s ) = 0.83 in the symmetry plane and over 
all space ||(a>||oo — 1.05. There is no anomalous negative vorticity and numerical solutions require less resolution. 



this is sucked underneath the primary vortex at t = 6. 
The t = 6 frame represents the vortices simulated with 
the 36th-order hyperviscosity [8, 11]. For t > 6 this sec- 
ondary vortex is dissipated by the hyperviscosity and cir- 
culation is dissipated, meaning these calculations are not 
faithfully representing the Euler equations. Without the 
hyperviscosity, numerical noise would dominate as an ex- 
tra boundary layer needs to be resolved. 

C. Final step 3) for purely positive 

What is apparently missing from the previous descrip- 
tion [3] is the addition of a mean shear designed to remove 
the final negative regions in the symmetry plane. This 
was achieved before [3] as part of an interpolation pro- 
cedure from a uniform mesh to a Chebyshev mesh. Here 
it is imposed. Details will appear in a full paper. Ini- 
tial vorticity in the symmetry plane and a slightly later 
time (t = 4.38) are shown. This is the initial condition 
for which we have now done up to 1024 x 256 x 2048 
calculations to assess the scaling proposed earlier [3]. 

V. A NEW APPROACH FOR DETERMINING 
WHETHER THERE IS SINGULAR BEHAVIOR 
OF THE PRIMARY PROPERTIES AND THE 
ASSOCIATED SCALING. 

Once reliable data (according to the criteria discussed 
above) has been obtained, it is common to interpret it 



in terms of power laws and other simple formulae. For 
example, assume that 

f(t) ~ C/{T C - ty . (4) 

To properly find all three free parameters (C, T c and 7) 
to a set of points requires a minimization procedure. 

Kerr (1993) avoided this by assuming particularly sim- 
ple values for 7 for several quantities. In particular 
7=1 was assumed for H^Hoo , for the maximum of the 
stretching of the vorticity (2) in the symmetry plane: 
max(a)|j /= o, and for the enstrophy production (3). This 
procedure was extended to the velocity by assuming that 
lu = 1/2 for sup(|u|) [5]. 

While fits with these assumptions gave consistent re- 
sults for the singular time T c , this consistency existed 
only at late times when resolution was becoming ques- 
tionable. Analysis of this data by two new methods has 
shown that the lack of scaling at earlier times is due in 
part to some of the more restrictive assumptions that 
were made. 



A. Three-parameter fitting 

Our first indication that earlier assumptions [3, 5] 
might be incorrect was obtained by allowing 7 to be free. 
The three parameters (C, T c and 7) were then obtained 
as follows: by minimizing the sum of squares of the dif- 
ferences between the logarithm of the data and the loga- 
rithm of the fit function, with respect to C and 7, allows 
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one to solve for these two parameters in terms of T c . Then 
the sum of squares is further minimized with respect to 
T c to obtain all three parameters. 

This analysis was applied to \\u>\\oo and Cl p , the enstro- 
phy production, both of which previously were assumed 
to have 1 / (T c — t) behavior. 

It was immediately observed that 

• The fitting parameters depended upon the range of 
times chosen. 

• 7 in each case was consistently greater than 1. 

This would not be inconsistent with known bounds. 
Recall if power law behavior is expected for ||u>||oo , (1) 
only requires that 7 > 1. 

B. Logarithmic time derivatives of enstrophy and 
||w||oo : instantaneous two-parameter fitting 

The original analysis would be possible if there is a 
secondary quantity which must go as 1/(T C — t) if the 
primary quantity obeys the power law 1/(T C — i) 1 . An 
example of a secondary quantity of that sort is the log- 
arithmic time derivative of the primary quantity, which 
can be computed if we know independently the quantity 
f(t) and its time derivative f(t). 

Therefore we propose a new approach to inferring a 
singular time and identifying the scaling behaviour: 

• Find a quantity f(t) whose growth and the growth 
of its time-derivative can be determined directly. 
Consider the new function: 

ff (t)=(|log/(t)) 1 = t= 1 -{T c -t). (5) 

Note that the parameter C drops out and the func- 
tion is linear, so we can predict instantaneous val- 
ues of 7 and T c by fitting this new function using 
adjacent points in time. 

• Calculating in this manner, using nearest neigh- 
bours in time, yields instantaneous predicted sin- 
gular times T c (t) and power laws j(t). These in- 
stantaneous parameters will generically depend on 
time. 

• See if these converge or relax (as time increases) . 

Pairs of quantities to which this procedure can be ap- 
plied are: 

• 1 1 a; I loo and its logarithmic time derivative a x (local 
stretching at the point Xqo). 

• Enstrophy ft and its production Q p (3) where we 
assume that 

Q ~ , Ci \ . (6) 



• Hclicity in the simulated quadrant of space and its 
production. 

This approach is applied to enstrophy and its pro- 
duction on our highest resolution simulations. Because 
enstrophy and its production are global quantities they 
converge numerically longer (to t = 11.25) than ||u>||oo > 
for which is less than 6 mesh points from the divid- 
ing plane for t > 10. Running estimates for T c and 70 
are shown in Fig. 2. For the new data (upper frames), 
the latest data point gives an estimate T c w 13.16 and 
70 ~ 0.47. The bottom frames show the same analy- 
sis applied to the data used by Kerr (1993) [3]. Nearly 
identical power laws are obtained (70 w 0.50), with a 
predicted singular time T c « 19.69, greater than in Kerr 
(1993). 

One advantage of finding running estimates of T c is 
that it can be used to identify cases that are not sin- 
gular, or would take an unusually long time to become 
singular. This is done by looking at the instantaneous 
estimated value of T c . If T c continues to increase with 
time, then there is evidence for regularity. In both the 
new calculations and for the data from Kerr (1993), even- 
tually the estimated T c decreases and relaxes to a finite 
value. It is quite possible that there is a large pre-factor 
in front of the power law, which the time dependence of 
the estimated 7n and Cq might be able to shed light on. 

This approach assumes smooth values for both quanti- 
ties in a pair. Unfortunately, we have found that because 
||w||oo sits on a steep gradient of a, values of on the 
lower resolution mesh were not smooth enough in time 
to perform this analysis. The analysis will be attempted 
on the higher resolution (||u>||oo ,&oo) data when that ad- 
ditional analysis of the new data sets is available. 

C. Convergence studies: Is the evidence for 
singularity conclusive? 

Further tests at higher resolution are needed to support 
the singular trends seen here. Both current cases (old 
and new data) could be reliably integrated up to times 
t w T c — 2.75. This would only be the beginning of the 
asymptotic regime of the potentially singular solution, as 
is suggested by the late-time behavior of the curves for 
the predicted singular time T c in Fig. 2. New calculations 
in progress should go beyond that barrier and help test 
the validity of the hypothesis of finite-time singularity. 

In this subsection we show resolution studies with n z 
fixed to give a flavor of what will be shown in the next 
paper (we have also made resolution checks with fixed 
n x or n y and varying the other two, not shown here). 
The resolution study is a classical tool to validate and 
find reliability times for the numerical results. Another 
now widely accepted study that we present is a spectral 
convergence test (Sulem et al., 1983 [12]) used recently 
by Cichowlas and Brachet (2004) (see [13] and references 
therein) , where the exponential decay of the energy spec- 
trum as a function of the wavenumber is employed to give 
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FIG. 2: Upper frames: Resolution study of the predicted singular time T c and the predicted exponent 7n in the power-law 
behavior of the total enstrophy using the new data at resolutions: 512 x 64 x 2048 (dashed), 512 x 128 x 2048 (dotted) and 
1024 x 256 x 2048 (solid). 

Lower frames: Predicted T c and 70. for the Kerr (1993) data at the highest resolution (solid). Dashed lines denote gaps in data. 
In the graphs for the predicted T c , the dash-dotted diagonal lines denote the T c = t singularity asymptote. 



a criteria for the reliability time. Finally we complement 
the above classical tests with our newly proposed tests of 
reliability: conservation of circulation through the sym- 
metry plane and conservation of circulation through the 
dividing plane. 

We consider first the behavior of local quantities, in 
the sense of Section III. Figure 3 (top) is a resolution 
study of the time dependence of the maximum vorticity; 
the bottom figure is a t = 10 anisotropic energy spec- 
trum E(k Xl t), defined by averaging the Fourier transform 
u(k',t) of the velocity field on flat duplicated sheets of 
width Ak x = 1, 

E{k x ,t) = \ I<i(k',t)| 2 . 

/c x -Afe x /2<|/c^|<fc ;c +Afc [c /2 

Following [13], we fit: log E(k x ) = C — nlog(k x ) — 28 k x 
(solid line in the figure). The test consists in monitoring 
the parameter 8 as a function of time. The idea is that 
(5, being the width in the complex plane of the analyt- 
icity strip of the velocity field, should always be numer- 
ically resolved, at least by the mesh size. Another way 



to look at this condition is to ask that the contribution 
of the exponential term to the change of log E(k x ) from 
the largest to the smallest scale allowed by the numer- 
ical resolution, be greater than a prescribed factor. In 
more explicit terms, we can only fully trust the simula- 
tion up until the condition 8 fc" lax > 1 is violated, where 
fc™ ax = n x /3 is the maximum relevant wavenumber of the 
Fourier representation. Notice that different authors use 
different factors in the RHS of the last inequality. For our 
t = 10 spectrum in resolution 1024 x 256 x 2048 we obtain 
(5fc™ ax w 1.07, and therefore our simulation is validated 
by this method up to t = 10. In this way we could extrap- 
olate the convergence of || us up to t = 10, whereas a 
conservative extrapolation based solely on the resolution 
study in figure 3 (top) would see the 1024 x 256 x 2048 
computation converged up to t = 9. 

We consider now the behavior of 2D integral quantities. 
Due to its 2D character, enstrophy in the symmetry plane 
flsp = J y=Q ujy dx dz, shown in figure 4 (top), is a more 
sensitive measure than total enstrophy £1 (eq.(3), figure 
not shown), which converges more rapidly than ilsp. A 
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FIG. 3: Top: Resolution study of jjtjjloo, for resolutions 
n x xn y x n z of: 512 x 64 x 2048 (dashed), 512 x 128 x 2048 
(dotted) and 1024 x 256 x 2048 (solid). Bottom: Anisotropic 
energy spectrum (direction k x ) at time t = 10 for resolution 
1024 x 256 x 2048. Points correspond to numerical data. The 
solid curve corresponds to the fit of the spectrum according 
to \ogE(k x ) = C — n\og(k x ) — 25 k x , where the fit interval is 
defined by the vertical dashed lines. 



conservative extrapolation would imply convergence of 
fl,sp up to times t J$ 11. 

Figure 4 (bottom) is a resolution study of the normal- 
ized error in the conservation of circulation through the 
dividing plane a z . We observe that, for a given resolu- 
tion, the numerically induced deviation in <j z becomes 
unstable after a certain time. Errors (and fluctuations 
thereof) less than 10~ 4 are acceptable, as long as they 
are stable. Then, a reasonable reliability time can be de- 
fined for each resolution as the time when the error in 
<7 2 attains its last extremum before the instability takes 
over. According to this criteria we conclude that the 
simulation at resolution 512 x 128 x 2048 (dotted line) 
is converged up to t « 10.7 and the simulation at res- 
olution 1024 x 256 x 2048 (solid line) converges up to 
t w 11.25. In order to display the unstable behavior of 
the mid-resolution simulation (dotted line) , we show data 
beyond its reliability time. 

Finally we return to Fig. 2, considering all three reso- 



FIG. 4: Resolution study of: the enstrophy in the symme- 
try plane Qsp (top); the error in the circulation through 
the dividing plane <r 2 normalized with the initial circula- 
tion through the symmetry plane a y (bottom), for resolu- 
tions 512 x 64 x 2048 (dashed), 512 x 128 x 2048 (dotted) and 
1024 x 256 x 2048 (solid). 



lutions. For each resolution, T c has a peak for t ~ 9—10, 
and then asymptotes, well within the reliability time for 
the highest resolution. Similar trends towards conver- 
gence appear for jq. The key question regarding the ex- 
istence of a finite-time singularity is if the curve for pre- 
dicted T c crosses the asymptote T c = t (dashed-dotted 
line) in a finite time or not, but we need further tests 
at higher resolutions (to be shown in a future paper) in 
order to conclude on these matters. 

We have enough resolution to conclude that the power 
laws are not the ones proposed in Kerr (1993,2005) [3, 5] 
but not enough resolution to reach definitive conclusions 
on the singular behavior, since ||w||oo does not converge 
as rapidly as the volumetric quantities studied. 



VI. GRAPHICS 

In this section, 3D isosurface contours of the vorticity 
modulus are shown, corresponding to the simulation of 
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FIG. 5: Euler anti-parallel vortices in full periodic domain near t = 2.51. Bright (yellow online) tubes are isosurface contours 
of vorticity modulus corresponding to 60% of the instantaneous maximum of vorticity modulus. Dark (red online) elongated 
blobs are isosurfaces corresponding to 90% of the maximum of vorticity modulus. 



initially anti-parallel vortices, using a resolution of 512 x 
128 x 2048 in the fundamental quarter of the full domain, 
corresponding to an effective resolution of 512 x 256 x 4096 
in the full domain. For memory-optimizing purposes, 
the output data used to make the figures has an effective 
resolution of 512 x 128 x 1024, corresponding to a memory 
size of 65 MB. The freeware visualization program Visit 
has been used to make the plots. 

Fig. 5 shows the vortices after some time of evolution 
in the whole periodic domain. The large tubes are iso- 
surfaces corresponding to 60% of the maximum vorticity 
modulus. These tubes cross along the periodic y-Axis 
and deform notably near the symmetry plane (y = 0). 
The elongated blobs in the interior of the tubes are iso- 



surfaces corresponding to 90% of the maximum vorticity. 
These isosurfaces are very localized and flattened near 
the symmetry plane. 

Fig. 6 shows successive snapshots at later stages of 
the flow, of isosurfaces corresponding to the following 
percentages of the instantaneous value of the maximum 
vorticity modulus: from outer to inner contours, 40%, 
60%, 80% and 90%. Only half of the total domain is 
shown so that a section of the isosurfaces through the 
symmetry plane (y — 0) is visible. The snapshots are all 
seen from the same angle and with the same zoom with 
respect to the fixed box. To read the snapshots going 
forward in time one advances from left to right and from 
top to bottom. 
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The flattening in the z-direction results in structures 
similar to flattened pillows with some curvature in x (and 
less in y) that becomes more pronounced at the later 
times along the cut at the symmetry plane. These em- 
pirical observations provide further support for the choice 
of anisotropic resolution in the simulations. 

VII. LOOKING FORWARD 

The calculation reported here was the largest possi- 
ble on the Warwick SGI Altix with our code. In the 
near future we will have a cluster capable of simulating a 
2048 x 1024 x 4096 mesh. We might also use UK national 
computing resources. 

Once our new cluster arrives we anticipate the follow- 
ing calculations: 

• Further tests to determine if an initial condition 
closer to that of Kerr (1993) [3] can be obtained. 

• After further resolution checks, at least one calcu- 
lation on a 2048 x 512 x 4096 mesh either on the 
new profile here or that of Kerr (1993) [3]. 

• At least one modest resolution calculation on the 
square-off initial condition of Fig. 1 . 

• Our goal in high resolution calculations will be to 
include spectral convergence tests, in particular the 
analyticity strip method (Sulem et al. [12], see also 
[13, 14] and references therein) which gives inde- 
pendent evidence of singular/nonsingular behavior 
of the flow and allows one to extrapolate the con- 
vergence Of || UJ | |oo . 

• Convergence of H^Hoo and other local quantities 
should allow us to study regularity bounds from 



Constantin, Fefferman and Majda [4] and from 
Deng, Hou and Yu [15]. 

The finite-time singularity hypothesis of the three- 
dimensional, incompressible Euler equations, leads to 
conclusions that are in qualitative agreement with Kerr 
(1993) [3]. However, we have found that the previously 
proposed scaling laws and estimated singular time must 
be modified. 

One possible outcome which will require further inves- 
tigation is whether constant circulation is trapped within 
the collapsing region. If this is confirmed, the two length 
scaling parameterization proposed [5] cannot be correct. 

It is anticipated that in addition to these much higher 
resolution anti-parallel calculations, there will soon be 
new high-resolution calculations of the Kida-Pelz flow 
(anticipated in these proceedings by the contribution of 
Grafke et al. [16]) and the Taylor-Green calculations 
initiated by Brachet et al. (1983) [17] will soon be 
continued by Brachet. If these prove to be singular, and 
anti-parallel, it is possible that they will reproduce the 
scalings hinted at here. 
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